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Abstract. We present studies of the collapse of neutron stars that undergo a hadron-quark 
phase transition. A spherical Lagrangian hydrodynamic code has been written. As initial 
condition we take different neutron star configurations taking into account its density, energy 
density and pressure distribution. The phase transition is imposed at different evolution times. 
We have found that a significant amount of matter on the surface can be ejected while the 
remaining star rings in the fundamental and first pressure modes. 



1. Introduction 

Black holes [BHs) and neutron stars (NSs) are certainly two major potential sources of 
gravitational waves {GWs). Unlike BHs, whose gravitational waveforms are specified essentially 
by their masses and angular momenta, the characteristics of the gravitational emission from NSs 
depend on the properties of the nuclear matter. 

When nuclear matter is compressed to a sufficiently high density, it turns into uniform three- 
flavor (n, d, and s) strange quark matter (SQM), since it is expected that SQM may be more 
stable than nuclear matter. The deconfined quark matter appears when the density is so high 
that the nucleons are "touching" each other. At this point, when the number density of nucleons 
n ~ 0(1 fm~^), the quarks lose their correlation with individual nucleons and appear in a 
deconfined phase. Since the density required for this to happen is not much higher than nuclear 
matter density (0.16 fm~^), the dense cores of neutron stars are the most likely places where 
the phase transition to quark matter may occur astrophysically. It should be noted that strange 
quarks (in a confined phase) could already exist in neutron stars with a hyperon coreP^. 

In principle, the existence of a thin crust of normal matter is possible at the surface of a 
strange star. On the other hand, if SQM is metastable at zero pressure, so that it is relatively 
more stable than nuclear matter only because of the high pressure in the cores of neutron stars, 
then the final products of the phase transition would be hybrid stars, which consist of quark 
matter cores surrounded by normal matter outer parts. 

The above arguments were made in [2] and those authors used polytropic equation of state to 
describe the complex nuclear matter structure. We proceed similar calculations with a realistic 
equation of state based on a field theoretical model. We also make use of a spherical Lagrangian 
hydrodynamical code instead of the Eulerian model used in . 



2. Hydrodynamical Collapse 

The spherical Lagrangian hydrodynamic code divides the star into several concentric shells. This 
number of shells is the optimal combination for the variable time step and the number of shells 
that gives the shortest processing time. Our results show that 1600 shells are enough to model 
satisfactorily all the stars here studied. To verify if the number of shells is an optimal choice 
we ran a couple of models with 6400 shells, and we see that the results are the same as for the 
models with 1600 shells. The shocks that appear are treated with the inclusion of the artificial 
viscosity of Richmeyer and Morton . 

The hydrodynamic equations that describe the dynamics of neutron stars are: 
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the mass conservation equation written in spherical coordinates; where p is the density, r the 
radial coordinate and v the velocity. In our code the continuity equation is calculated directly 
from the grid. The specific volumes of the Lagrangian shells give the density of the star. 
The equation of motion is given by: 
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where the velocity is v; the pressure of the nuclear matter is given by the model described in 
the next section. The gravitational potential is (p. The energy equation is written as: 

DE _PDp 
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where E is the thermal energy per gram. 

3. Neutron Star vs. Equation of State 

In our investigation we consider a phenomenological QHD Lagrangian density [5l [6l [7] which 
contains the fundamental baryon octet, the scalar meson fields a, /o and the vector meson fields 
w, /9, (j). Additionally, it contains the lightest charged leptons (e and p) to allow for charge 
neutrality, 
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where the parameterized coupling constants [5] of the theory are defined as 

9^B = ^*BI3 9ujb 
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with 
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For the quark matter the MIT Lagrangian was adopted. Here we have assumed that 
the phase transition is of first order. This last possibihty is not favored energeticahy when the 
interface tension between quarks and hadrons is important [9l[l0]. Thus, the physical conditions 
at deconfinement were estimated from Gibbs criteria. 

4. Numerical Results and Conclusions 

We have obtained some important results which seems to be in complete agreement with those 
obtained by p]. The equation of state is shown in Fig.l for three different parameter a which 
leads to different hadron-quark phase transition pressure and density. Three stages of evolution 
are shown in Fig. 2. The star radius shrinks and the central density is increasing during the 
collapse. In Fig. 3 one can find the density evolution of each shell. The oscillation of the inner 
shells is clear in the final stage of evolution as it is the ejection of the external shells. The fourier 
transform is shown in Fig. 4. This result shows a peak very close to the one shown by [2], at 
3.1kHz, which shall correspond to the fundamental mode as well as other peaks, corresponding 
to the first pressure modes. It is also easy to find that most of the energy is released in the 
fundamental mode, as always claimed. In Fig.5 the radius and density evolution are plotted 
against time in this same figure. 

Finally, after a future detection, one can try to fit the nuclear matter parameters to obtain the 
frequency of the detection, remembering that, if the frequency falls on the 3. 0-3. 4kHz bandwidth, 
the Brazilian Schenberg antenna will be able to identify the direction of the source, as shown in 
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Figure 3. Density shells evolution from collapse to oscillation modes. Time is expressed in 
seconds. 
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Figure 4. Fourier transform of the oscillations. Frequency is written in kHz. 



